Can high energy neutrino annihilation on relic 
neutrinos generate the observed highest 
energy cosmic-rays? 



E. Waxman 

Institute for Advanced Study, Princeton, NJ 08540; e-mail: waxman@sns.ias.edu 



Abstract 

Annihilation of high energy, ~ 10 21 eV, neutrinos on big bang relic neutrinos of 
~ leV mass, clustered in the Galactic halo or in a nearby galaxy cluster halo, 
has been suggested to generate, through hadronic Z decay, high energy nucleons 
and photons which may account for the detected flux of > 10 20 eV cosmic-rays. We 
show that the flux of high energy nucleons and photons produced by this process is 
dominated by annihilation on the uniform, non-clustered, neutrino background, and 
that the energy generation rate of ~ 10 21 eV neutrinos required to account for the 
detected flux of > 10 20 eV particles is > 10 48 erg/Mpcyr. This energy generation 
rate, comparable to the total luminosity of the universe, is ~ 4 orders of magnitude 
larger than the rate of production of high energy nucleons required to account for the 
flux of > 10 19 eV cosmic-rays. Thus, in order for neutrino annihilation to contribute 
significantly to the detected flux of > 10 20 eV cosmic-rays, the existence of a new 
class of high-energy neutrino sources, likely unrelated to the sources of > 10 19 eV 
cosmic-rays, must be invoked. 
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1 Introduction 



The Fly's Eye [1] and AGASA [2] experiments confirmed the existence of a 
break in the energy spectrum of high energy cosmic rays at ~ 5 x 10 18 eV, 
for which evidence existed with weaker statistics in the data of previous ex- 
periments (Haverah Park, Yakutsk, Sugar, see e.g. [3] for a review). Fly's Eye 
data also strengthen the evidence for a change in primary composition from 
predominantly heavy nuclei below the break to predominantly light nuclei 
above the break. These features strongly suggest, when coupled with the lack 
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of anisotropy that would be expected for cosmic-rays (CRs) of Galactic origin, 
that below ~ 10 19 eV the CRs are mostly heavy ions of Galactic origin, and 
that an extra-Galactic component of protons dominates above ~ 10 19 eV. This 
conclusion is further supported by the fact that the CR energy spectrum is 
consistent with a cosmological distribution of sources of protons, with injec- 
tion spectrum dN/dE oc E~ 2 ' 2 typically expected for Fermi acceleration [4]. 
In particular, there is evidence for the existence of a Greisen-Zatsepin-Kuzmin 
(GZK) "cutoff", i.e. for the suppression of CR flux above ~ 5 x 10 19 eV ex- 
pected due to interaction of protons with the microwave background radiation 
[5]. 

The evidence for GZK suppression is strengthened by recent AGASA data [6]. 
In Fig. 1, the CR spectrum reported by the Fly's Eye and the AGASA experi- 
ments [1,6] is compared with the flux expected for a homogeneous cosmological 
distribution of sources, each generating a power law differential spectrum of 
high energy protons dN/dE oc E~ 2 ' 2 (For the model calculation we have used 
a flat universe with zero cosmological constant, Hq = 75km s _1 , and time inde- 
pendent energy generation rate per comoving unit volume 5 x 10 44 erg/Mpc 3 yr; 
The spectrum is insensitive to the cosmological parameters and to source evo- 
lution, since most of the cosmic rays arrive from distances < lGpc [4]). The 
deficit in the number of events detected above 5 x 10 19 eV, compared to a 
power-law extrapolation of the flux at lower energy, is consistent with that 
expected due to a cosmological GZK suppression. However, with current data 
the "cutoff" is detected with only 2a significance [7]. 

The number of events detected above 10 20 eV is consistent with that expected 
based on the cosmological model presented in Fig. 1 (There is an apparent 
"gap" between the highest and second highest energy events detected by the 
Fly's Eye [8]. However, assuming that the cosmological model is valid, the 
probability that such an apparent "gap" would be observed is ~ 15% [4,8]). 
Nevertheless, the detection of > 10 20 eV events does pose challenges to most 
models of CR production. The high energies rule out most of the accelera- 
tion mechanisms so far discussed [9], and since the distance traveled by such 
particles must be smaller than lOOMpc [10] due to their interaction with the 
micro-wave background, their arrival directions are inconsistent with the po- 
sition of astrophysical objects, e.g. jets of powerful radio galaxies [11], that 
are likely to produce high energy particles [12]. 

Cosmological 7-ray bursts (GRBs) are likely sources of high-energy CRs, which 
may account for the CR flux above 10 19 eV as well as for the > 10 20 eV events 
[13]. This model recently gained support form GRB afterglow observations 
[14] . Other models for the production of ultra-high energy CRs were suggested, 
where the highest energy events are produced by the decay of super-massive 
elementary particles related to grand unified theories (see, e.g., [15] for recent 
review) . Sources of such particles may be topological defects, left over from a 
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phase transition associated with the symmetry breaking of the grand unified 
theory [16]. While no firm prediction exists of the CR flux in these theories, a 
generic feature of the super-massive particle decay scenarios is that the injec- 
tion spectrum is much harder than expected for Fermi acceleration. Therefore, 
this scenario can account only for the flux of > 10 20 eV particles, and can not 
simultaneously explain the origin of 10 19 -10 20 eV CRs. 

It has recently been suggested that annihilation of high energy, ~ 10 21 eV, 
neutrinos on big bang relic neutrinos of ~ leV mass, clustered in the Galactic 
halo or in a nearby galaxy cluster halo, may generate high energy nucleons and 
photons which may account for the detected flux of > 10 20 eV cosmic-rays [17]. 
The existence of > 10 21 eV neutrino flux was argued plausible based on the 
argument that the mechanism producing the observed high energy, > 10 19 eV, 
particles, most likely protons, also produces charged pions of comparable en- 
ergy, which subsequently decay to produce neutrinos. It was suggested that 
the generation spectrum extends well beyond 10 20 eV, and that while nucleons 
produced by a distant source, e.g. a powerful radio galaxy, lose their energy 
interacting with the micro-wave background, high-energy neutrinos propa- 
gate without energy losses and may annihilate on relic neutrinos, producing 
> 10 20 eV nucleons and photons at small distances that would allow them to 
propagate to Earth. In Sec. 2 we derive the energy density of high energy 
neutrinos required to account for the observed rate of > 10 20 eV air showers. 
The implications of our results are discussed in Sec. 3. 



2 The high energy neutrino density 

Hot Big Bang cosmology predicts the existence of a neutrino background, 
similar to the photon micro- wave background. The energy distribution of neu- 
trinos (of mass m v below lMeV, the decoupling temperature) is predicted to 
be a Fermi-Dirac distribution with temperature smaller than that of the pho- 
ton background, T 7 = 2.73K, by a factor (11/4) 1 / 3 . The corresponding number 
density of neutrinos is % = 54 cm~ 3 , with a similar density of anti-neutrinos 
(we assumed a single helicity state; the presence of two helicity states would 
increase the density by at most a factor 2, and will not change the conclu- 
sions derived below). Neutrino mass smaller than ~ 5eV is consistent with Big 
Bang cosmology [18]. Furthermore, there are several indications that neutri- 
nos are indeed massive. The simplest explanation for the discrepancy between 
observed and predicted solar neutrino fluxes (see [19] for a review) involves 
u e — oscillations driven by neutrino mass difference, and the common ex- 
planation to the atmospheric neutrino anomaly involves — v T oscillations. 
While the mass difference implied for v e — oscillations is much smaller than 
leV 2 , the common explanation to the atmospheric anomaly is oscillation to 
v T with mass ~ 0.03eV [20]. 
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High energy neutrinos may interact with the background neutrinos and anni- 
hilate, producing Z-bosons which immediately decay with 70% probability for 
hadronic decay. With standard weak interaction, the annihilation cross sec- 
tion is strongly peaked at the resonant neutrino energy e R = M\j2m v = 
4 x 10 21 (my/leV) _1 eV, and the energy-averaged cross section, defined as 
a = / dsa vv {s)jM\ where s is the square of the center of momentum en- 
ergy, is a = AttGf/V^ = 4 x 10~ 32 cm [21]. The relative energy width of the 
resonance peak (FWHM) is ~ 3%. Z decay produces on average 2.7 nucleons 
and anti-nucleons, and 30 high energy photons through 7r° decay. A single nu- 
cleon carries, on average, a fraction ~ 0.025 of the Z energy [22], and a single 
photon carries, on average, energy lower by a factor of ~ 10 compared to the 
nucleon energy [17,22]. For a neutrino mass m u < lOeV for the background 
neutrinos, the highest energy nucleons produced by the decay of Z-bosons cre- 
ated from annihilation of resonant (e R ) high energy neutrinos and background 
neutrinos have energy > 2 x 10 20 eV, and could in principle account for the 
observed highest energy CR events. 

Let us first consider the rate of high energy showers due to annihilation of 
high energy neutrinos with a homogeneous neutrino background. The rate per 
unit volume of Z production is 



where % is the number density of background neutrinos (or anti-neutrinos) 
and dn/de is the number density per unit energy of high energy neutrinos. We 
note here that, as suggested by oscillation evidence, neutrino flavor eigenstates 
may be different than their mass eigenstates, and since neutrinos produced by 
astrophysical sources are expected to be electron and muon neutrinos from 
pion decay, the Z production rate is smaller than given by (1) by a factor 
I < v a \v m > | 2 , where |z/ a > is the flavor eigenstate and \u rn > is the massive 
neutrino eigenstate. If, as suggested by present data, the massive, ~ leV, 
neutrino is mainly mixed with u T , then the annihilation rate of v e ^ expected 
to be produced by astrophysical sources will be smaller than given in (1). 
However, since we are interested in obtaining a lower limit to the density of 
high energy neutrinos required to produce hz high enough to account for the 
observed flux of CRs above 10 20 , we conservatively assume that | < v a \v m > | 2 
is of order unity for the high energy neutrinos. 

The integral in (1) is dominated by the contribution from e ~ e R . We may 
therefore replace the derivative dn/de in the integrand with its value at e = e R . 
For m u ~ leV the background neutrinos are nonrelativistic, and the square 
of the center of momentum energy is s = 2em u c 2 . Changing the integration 
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variable e to s, we find 



nz = cubct 



dn 
de 



(2) 



The high-energy protons and photons produced by Z decay lose energy as 
they propagate through the microwave background. This limits the distance 
out to which Z decays contribute to the observed rate of air showers above 
10 20 eV to tcr < lOOMpc. The flux per solid angle of high energy cosmic-rays 
producing showers above 10 20 eV, and originating in Z decays, is therefore 



Jcr = -^n z N CR (m u )r C R, 



(3) 



where Ncii{m u ) is the average number of protons and photons above 10 20 eV 
produced by the decay of a Z-boson of energy e R . 

Using eqs. (2,3) we may derive the number density [edn/ de] e=e R of high energy 
neutrinos in the energy range Ae ~ around e = e R , required to produce a 
given CR flux j CR , 



e=e R cnB<?rcRNcR{m y ) ' 



n 



R 



dn 
de 



The total energy density associated with high energy neutrinos is expected 
to be larger by a factor j\ e > 1 than the energy e R n R associated with neu- 
trinos in the energy range Ae ~ e R , since one does not expect the neutrino 
energy distribution to be strongly peaked at e R . If the high energy neutrino 
flux is indeed produced by the sources of > 10 19 eV CRs, then one expects 
the neutrino distribution to be approximately given by dn/de oc e~ 2 [23], 
for which the energy density in neutrinos of energy 10 19 eV < e < is 
E u (> 10 19 eV) = e R n R log(e R /10 19 eV) ~ 5e R n R . Thus, the energy density 
of > 10 19 eV neutrinos required to produce the observed rate of showers above 
10 20 eV, j CR ^ 10- 16 m- 2 s~ 1 sr- 1 (see Fig. 1), is 



E v (> 10 iy eV) = f Af 



AncM 2 z j C R 



n B o-rcRm u N C R(m u ) 



* 10- lb (W5) 



N CR {m v )m v tcr 
3eV lOOMpc 



/ 3 

erg/ cm 



(5) 



E u depends only weakly on m u , since N C R(m u )m u ~ 3eV independent of the 
value of m u . For m u ~ leV, the resonant energy is e fi ~ 4 x 10 21 eV, and 
the average energy of the nucleons produced in the Z decay is ~ e R /A0 ~ 
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10 20 eV. Since the average number of nucleons produce in the decay is 2.7, we 
have NcR(m u )m u ~ 3eV for m u ~ leV. For higher m v the resonant energy 
decreases, and only a small fraction of the resulting nucleons have energy above 
10 20 eV. For m u ~ lOeV, e fi ~ 4 x 10 20 eV and since (on average) only ~ 10% 
of the nucleons carry ~ 1/4 the Z energy [22], we have NcRim^m,, ~ 3eV for 
m u ~ lOeV. For m v ~ O.leV, the resonant energy is e R ~ 4 x 10 22 eV, so that 
the average energy of the photons produced in the decay is ~ 6^/400 ~ 10 20 eV, 
and the photons also contribute to jcn- Since the average number of photons 
produced in the decay is 30, we have N C R{ra u )m u ~ 3eV for m v ~ O.leV. 

The energy density E V) required for Z decays to contribute significantly to the 
observed rate of > 10 20 eV air-showers, may be reduced if the local density of 
background neutrinos is higher than the average density tir, due to cluster- 
ing. For m u ~ leV, the background neutrinos are non-relativistic at present, 
and are expected to cluster in the gravitational potential wells of galaxies and 
clusters. Let us first consider clustering in the halo of the Galaxy. The density 
of neutrinos in the Galactic halo depends on the details of the halo forma- 
tion. However, Pauli's exclusion principle allows to put an upper limit to the 
number density ric of neutrinos clustered in the Galactic halo. This limit is 
approximately given by n c < (47r/3)(p max //i) 3 , where p max ~ 2 l l 2 m u v r is the 
maximum momentum of bound neutrinos and v r = 220km/s is the Galactic 
rotation speed. A limit lower by a factor of 2 is obtained by requiring the 
neutrino phase space density, approximately given by nc/ (47rp 3 iax /3), to be 
smaller than the maximum phase space density of the uniform background, 
1/2/i 3 . This requirement follows from the assumption that neutrino cluster- 
ing results from gravitational collapse of the uniform neutrino background. In 
this process, the phase space density is conserved along particle trajectories, 
and therefore the maximum phase space density can not increase [24]. The 
neutrino density in the halo is therefore constrained by 

n c < 1.5 x 10 3 ( ^-Y ( V ; , V cm" 3 . (6) 
° VleVy V220km/s/ V ; 

Using eqs. (2) and (3), replacing hr with nc and tcr with the halo radius 
tr, we find that the ratio of the flux of > 10 20 nucleons and photons due to 
Z production by annihilation with neutrinos clustered in the Galactic halo to 
the flux due to annihilation on the uniform background is 

/G//B<0 .0l(^) 3 ^. (7) 

Thus, for m u ~ leV, the rate of > 10 20 eV showers due to annihilation on 
neutrinos clustered in the Galactic halo is two orders of magnitude smaller 
than that due to annihilation on the uniform neutrino background. The con- 
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tribution of Galactic halo neutrino annihilation increases with m u , and for 
m v = lOeV may be larger than the contribution due to annihilation on the 
uniform background by a factor ~ 10. 



Finally, let us consider the contribution to jcn due to annihilation on back- 
ground neutrinos clustered in a nearby galaxy cluster. The rate of Z produc- 
tion in the cluster is 



N Z = Sj^ac 



dn 
de 



(8) 



where M c is the cluster mass, and f v is the fraction of cluster mass contributed 
by (background) neutrinos. The cluster mass is given by its velocity dispersion 
a v and radius R, M c = 2Ra 2 / G. The ratio of high energy shower rate due to 
annihilation in the cluster, NzNqr/ '^Tid 2 where d is the cluster distance, to 
that due to annihilation on the uniform background, nz^cBj'CR-, is 

fc/fB * °-° 3 (leV) (20^) b3l^(o^bInVi) • (9) 



Clearly, the contribution to high energy showers from neutrino annihilation in 
a nearby galaxy cluster is much smaller that the contribution from annihilation 
on the uniform background. We note here that in deriving eq. (8) we have 
assumed that the neutrinos are non-degenerate, so that the fraction f u is not 
limited by Pauli's principle. This is not valid for small m v . The cluster density 
at radius R is <j 2 /2iiGR 2 , implying a neutrino (anti-neutrino) density 

n c = 2, x 10U/0.1) (J^y (^y 2 «n- (10) 



This density is lower than the upper limit imposed by eq. (6) for 

m „ > 0.7fA) I/4 f_«_)- I/, (_^_V 1/4 .v, (11) 

V - 1 / \iMpc; v 600km / s / 



where we have used v r = 2 l / 2 a v . For m„ < 0.7eV the neutrino density is 
limited by the exclusion principle, and the neutrino fraction of cluster mass is 
proportional to m\. Thus, although for low m u , m„ ~ O.leV, photons produced 
in the Z decay may contribute to the particle flux above 10 20 eV, the decrease 
in f v implies that the ratio fc/fB is increasing with m u for m„ < 0.7eV. 
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3 Implications 



We have shown that the energy density of high energy neutrinos required to 
generate the observed flux of > 10 20 CRs by producing high energy nucleons 
and photons via the decay of Z-bosons resulting from the annihilation of the 
high energy neutrinos on a uniform neutrino background is ~ 10~ 16 erg/cm 3 , 
independent of neutrino mass m„ for O.leV < m v < lOeV (cf. eq. 5; the 
required energy density is higher for m v < O.leV, and for m v > lOeV the 
Z decay products are not energetic enough to account for the observed > 
2 x 10 20 eV CRs). The > 10 20 CR flux produced by annihilation on background 
neutrinos clustered in a nearby galaxy cluster is 2 orders of magnitude smaller 
than that due to annihilation on the uniform background, independent of 
the value of m v (cf. eq. (9)). For m v ~ leV, the CR flux resulting from 
annihilation on neutrinos clustered in the halo of our Galaxy is smaller by at 
least a factor ~ 100 than that due to annihilation on the uniform background 
(cf. eq. (7)). The relative contribution of Galactic halo neutrino annihilation 
increases with increasing m u , and may become comparable to the uniform 
neutrino background contribution for m v ~ 5eV. For m v ~ lOeV, annihilation 
on Galactic halo neutrinos may dominate, giving a CR flux higher by a factor 
~ 10 compared to the flux due to annihilation on the uniform background. 
In this case, the lower limit to the energy density of high energy neutrinos 
required to produce the observed flux of > 10 20 eV CRs is lowered by a factor 
of ~ 10 to ~ 10~ 17 erg/cm 3 . 

For m u ~ leV, the energy density in high energy neutrinos, ~ 10~ 16 erg/cm 3 , 
required to account for the observed flux of > 10 20 eV CRs is ~ 4 orders 
of magnitude higher than the energy density in observed > 10 19 eV cosmic- 
rays, A-kJE/c ~ 10 _20 erg/cm~ 3 (see Fig. 1). It is important to note here that 
although the life time of high energy neutrinos is longer than that of high 
energy protons, since 10 19 eV protons suffer energy loss due to pair production 
on the microwave background while neutrinos lose energy only due to redshift, 
the life time of a 10 19 eV proton, ~ 3 x 10 9 yr, is comparable to the Hubble time, 
~ 10 10 yr, and therefore to the neutrino life time. Therefore, the average (over 
time and volume) energy generation rate of high energy neutrinos required to 
produce the neutrino energy density necessary to account for the > 10 20 eV 
events is at least ~ 4 orders of magnitude higher than the generation rate 
~ 5 x 10 44 erg/Mpc 3 yr of > 10 19 eV protons required to account for the observed 
> 10 19 eV CR flux [4]. This implies that in order for neutrino annihilation 
to contribute significantly to the detected flux of > 10 20 eV cosmic-rays, the 
existence of a new class of high-energy neutrino sources, likely unrelated to the 
sources of > 10 19 eV cosmic-rays, must be invoked. Furthermore, the energy 
generation rate of the high energy neutrino sources must be ~ 10 49 erg/Mpc 3 yr, 
comparable to the total photon luminosity of the universe. 
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Fig. 1. The CR flux expected for a homogeneous cosmological distribution of 
sources, each generating a power law differential spectrum of high energy protons 
dN/dE oc E- 2 - 2 , compared to the Fly's Eye and AGASA data (The AG AS A flux 
at 3 — 10 x 10 18 eV is ~ 1.7 times higher than that reported by the Fly's Eye, cor- 
responding to a systematic ~ 20% larger estimate of event energies in the AGASA 
experiment compared to the Fly's Eye experiment [1,2]; We have therefore shifted 
upward (downward) the Fly's Eye (AGASA) event energies by 10%). Integers indi- 
cate the number of events observed, la energy error bars are shown for the highest 
energy events. The dashed line denotes the power law fit by Bird et al. [1] for the 
extra-galactic flux dominating above ~5x 10 18 eV. 
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